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We investigate numerically the zero-temperature physics of the one-dimensional Bose-Hubbard 
model in an incommensurate cosine potential, recently realized in experiments with cold bosons in 
optical superlattices [L. Fallani et al, Phys. Rev. Lett. 98, 130404, (2007)]. An incommensurate 
cosine potential has intermediate properties between a truly periodic and a fully random potential, 
displaying a characteristic length scale (the quasi-period) which is shown to set a finite lower bound 
to the excitation energy of the system at special incommensurate fillings. This leads to the emergence 
of gapped incommensurate band-insulator (IBI) phases along with gapless Bose-glass (BG) phases 
for strong quasi-periodic potential, both for hardcore and softcore bosons. Enriching the spatial 
features of the potential by the addition of a second incommensurate component appears to remove 
the IBI regions, stabilizing a continuous BG phase over an extended parameter range. Moreover we 
discuss the validity of the local-density approximation in presence of a parabolic trap, clarifying the 
notion of a local BG phase in a trapped system; we investigate the behavior of first- and second- 
order coherence upon increasing the strength of the quasi-periodic potential; and we discuss the 
ab-initio derivation of the Bose-Hubbard Hamiltonian with quasi-periodic potential starting from 
the microscopic Hamiltonian of bosons in an incommensurate superlattice. 
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Localization effects in random potentials represent a 
striking manifestation of the wavelike nature of quantum 
particles, and they are particularly pronounced in dimen- 
sions less than three, where, in absence of interaction, 
the suppression of transport due to disorder occurs at all 
energy scaleS^^^l, A particularly intriguing aspect of local- 
ization phenomena is their subtle interplay with particle- 
particle interaction. In the following we will specify our 
discussion to bosons. The screening effect of weak inter- 
actions can lead to a transition from an localized (Ander- 
son) insulator to a superfluid, namely t o th e persistence 
of superfluidity in a disordered potential^. Systems on 
a commensurately filled lattice without disorder exhibit 
a gapped Mott-insulating (MI) phase for strong interac- 
tions, competing with a superfluid (SF) phase as the in- 
teraction is reduced; in presence of disorder, the gapless 
Bose-glass (BG) insulator competes with the MI when 
the dis order strength becomes comparable with the Mott 
garPEI. 

Recent experiments on trapped cold atoms have 
demonstrated the ability of loading an a tomic Bose gas 
in a well controlled disordered potentiaPEHMMIIl 
particular, the regime of strong repulsion in a strong 
psewdo-disordered potential has been achieved by load- 
ing cold bosons in a deep one-dimensional quasi-periodic 
(QP) optical superlattice, formed by two standing waves 
with incommensurate wavelenghtgi^. In the insulating 
phase of the system, Bragg spectroscopy reveals an ex- 
citation spectrum which strongly differs from that of a 
MI, showing that the QP potential introduces many new 
states in the spectral Mott gaps, caused by the strong 
repulsion. This feature indirectly suggests that the in- 
sulating phase measured in the experiments could be a 
BG; yet, direct experimental inspection in the low-energy 
spectrum, which is decisive to distinguish a BG insulator 
from a more conventional MI, is currently not available. 



On the theoretical side, the study of tight-binding 
models in QP potentials has a long historyi^^'^ 
connected to the fundamental question of the fate of An- 
derson localization in quasi-random potentials. Although 
Bloch's theorem docs not apply to QP potentials, for one- 
dimensional systems it was soon realized that the single- 
particle spectrum of the system is composed of bands, 
made of either all extended or all localized states de- 
pending on the potential strength. Therefore, filling such 
a system with spinless fermions leads to the appearence 
of incommensurate band-insulator (IBI) phases alternat- 
ing with metallic phases for weak QP potential, and with 
Anderson-insulating ones for strong QP potential. Re- 
markably, the very same picture can be ex actly recovered 
for bosons in the hardcore repulsive limUP^EH^ with the 
following correspondence between fermionic and bosonic 
phases: "metallic" — > "superfluid" and "Anderson insu- 
lator" "Bose glass" . 

In this paper we provide a systematic study of one- 
dimensional bosons in a QP potential. We explicitly 
map out the phase diagram of the system in the hard- 
core limit through exact diagonalization, reconstructing 
the tight alternation of the SF, IBI and BG phases upon 
changing the chemical potential. Making use of quan- 
tum Monte Carlo we can then move on to the softcore 
case which is most appropriate for the experiments^. 
There we find that a strong QP potential, equalling in 
strength the inter-particle repulsion, completely removes 
the MI phase, leaving space for various insulating phases 
at fractional fillings: similarly to the hardcore case, slight 
changes of the model parameters (chemical potential and 
hopping amplitude) can drive the system from BG to IBI 
in a tight alternation. Therefore pseudo-disorder created 
by a single incommensurate potential component leads to 
a significantly different p hase dia gram than in the case of 
a truly random potentiaPSIMII]^ where a continuous BG 
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appears, becoming the only insulating phase for strong 
enough disorder. The fundamental differences between 
a truly random potential and a QP one are explained 
quantitatively within a generalized atomic-limit approxi- 
mation, which captures how the statistics of spatial struc- 
tures in the potential reflects itself in the spectral prop- 
erties. On the basis of this observation, we extend the 
degrees of freedom of the QP potential beyond what has 
been so far explored in experiments^, in search for the 
minimal setup which would lead to a similar behavior to 
the case of a truly random potential. Wc find that the ad- 
dition of a further incommensurate component, leading 
to a 2-color QP potential, causes the removal of the IBI 
phases, and hence to the appearence of a BG persisting 
over an extended filling range. 

Moreover we discuss the behavior of the system in pres- 
ence of a parabolic trap, mimicking the experimental sit- 
uation encountered in current optical-lattice setups. The 
local-density approximation (LDA), generally successful 
in the absence of a QP potential, is critically discussed 
in its presence. On the one hand the density profile and 
local particle-number fluctuations are correctly predicted 
by LDA; on the other hand the identification of the lo- 
cal behavior in the trapped system with that of the bulk 
system at the same chemical potential does not generally 
hold, due to the severeness of the finite-size effects im- 
posed by the trap in presence of the QP potential. At the 
same time we show that a 2-color QP potential can real- 
ize a uniform BG phase in the trap, on average over the 
fluctuations of the spatial phases of the potential com- 
ponents. Finally we make a closer connection to optical 
lattice experiments by investigating behavior of the cur- 
rently accessible observables, which are shown to reveal 
the full evolution of phase and density correlations in the 
system upon changing the strength of the QP potential, 
and we discuss in details the ab-initio derivation of the 
one-dimensional Bose-Hubbard Hamiltionian in the case 
of an incommensurate optical superlattice. 

This paper is structured as follows. Section |l] intro- 
duces the Bose-Hubbard model in an incommensurate 
cosine potential, and the numerical methods used in the 
paper; Section|ll]discusses the phase diagram in the hard- 
core limit; the phase diagram of the softcore model is 
then presented in Section Section |IV| discusses the 
emergence of the IBI in presence of the QP potential; 
Section |V] focuses on the local-density approximation in 
the trapped system; Section [VI| shows results for the ex- 
perimental observables; and Section [Vll| presents the ab- 
initio derivation of the Bose-Hubbard Hamiltonian. Sum- 
mary and conclusions are presented in Section [VIII| 



I. MODEL AND METHODS 
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where bi, b] are bosonic operators, J is the hopping am- 
plitude, U is the on-site repulsion, V2 is the strength of 
the QP potential, and fi the chemical potential. The 
incommensuration parameter < a < 1 is in principle 
an irrational number; making use of periodic boundary 
conditions, it is appropriate to take a = Ncydcs/L (with 
non-integer L/7Vcycios) so that the incommensurate cosine 
potential has a period exactly equal to the lattice size L. 

The realization of the Bose-Hubbard model by loading 
cold bosons in the lowest band of an optical lattice^ 
has been demonstrated by now in a large variety of 
experiments'^"^ ■^'^ . More recently Fallani et aZ.^'^ have 
been able to add a QP potential to a one-dimensional 
optical lattice by application of a secondary standing 
wave, with an incommensurate wavelength A2 with re- 
spect to that (Ai) of the primary lattice. If the inten- 
sity of the secondary wave is significantly weaker than 
that of the primary lattice, the main effect of the A2- 
lattice is to modulate the depth of the potential wells 
of the primary Ai-lattice, as modeled by the V2-term in 
Eq. ^ (see Section VII for a thorough discussion of 
this point). We choose the representation of the incom- 
mensurate cosine potential as in Eq. ^ such that the 
potential strength V2 reflects directly the intensity of a 
secondary standing wave at wavelength A2 — aXi: given 
that the first, dominant standing wave "discretizes" the 
space to points Xi = (Ai/2) i, the secondary standing 
wave creates a potential on such points which is propor- 
tional to cos^[(27r/Ai-|-(/))a;i] = cos^(a7ri-|-0) as in Eq. (|2|. 
The wavelength relation defines the incommensuration 
parameter, a = A1/A2 = 830.7/1076.8 = 0.77137... . In 
the following we will use different rational approximants 
thereof {a « 97/126, 830/1076). 

A fundamental limit of the Bose-Hubbard model is the 
hardcore case C/ — s- 00, in which double occupancy is 
strictly forbidden, and which has been experimentally 
demonstrated in Ref. 26. In the hardcore limit the con- 
straint of forbidden double occupancy can be directly in- 
corporated in the operator algebra by introducing on-site 
anticommutation rules for the hardcore-boson operators, 
{ai,al} = 1, {a-^'', a|^'} = 0, coexisting with ordinary 
bosonic commutation rules between different sites. Hard- 
core bosons can be exactly mapped onto spinless fermions 
fi, fj through the Jordan- Wigner transformatiorP3 



We investigate the one-dimensional Bose-Hubbard 
model in an incommensurate cosine potential, with 
Hamiltonian 
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so that the resulting Hamiltonian is that of free 
fermions in a potential: 



scope of the present paper, and we postpone it to a forth- 
coming publication. 
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The boundary phase term resulting from the non-local 
nature of the Jordan- Wigner transformation reads 9 — 
[{N -I- 1) ± IJtt, where N is the particle number and the 
sign -I- (— ) applies to periodic (antiperiodic) boundary 
conditions. 

Once reduced to spinless fermions, the hardcore bo- 
son problem is exactly solvable through simple single- 
particle diagonalization. In the following we will focus 
our attention on two fundamental quantities, namely the 
superfluid density, estimated through the energy differ- 
ence between periodic (4-) and antiperiodic (— ) boundary 
conditions. 



L 



(5) 



and the compressibility k = dn/dji, where n = N/L is 
the particle density. 

In the softcore case of the original Hamiltonian Eq. ([T]) 
an exact solution is not available. In this case we 
make use of quantum Monte Carlo (QMC), based on 
the Stochastic Series Expansion formulation and on the 
directed-loop algorithrrPS The simulations are system- 
atically done at temperatures low enough (typically T ~ 
J/L) so as to remove significant thermal effects. The 
truncation of the local Hilbert space up to a maximum 
bosonic occupancy n„iax is accurately done to avoid sig- 
nificant truncation errors (we went up to rimax = 6 for 
the largest fillings considered, n ^ 3). When averag- 
ing over different realizations of the QP potential, we 
have considered typical samples of 100 — 300 realiza- 
tions. The QMC method we use is formulated in the 
Fock-state basis, so that all quantities which are diag- 
onal in that basis (density, compressibility, etc.) are 
straightforwardly evaluated. The superfluid density is es- 
timated through the fluctuations of the winding number 
W in the worldline configurations produced during the 
simulatioiP^, ps = (VF^)L/(2/3 J). Finally we calculate 
the momentum distribution 



{h]hra) 
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obtained through the statistics of the directed-loop 
updatePI. Most of the QMC results of the paper are 
based on a standard grand- cano nical algorithm, alre ady 
well documented in the literatur^2ll30l_ Sections |v] |v| 
and VII we also present results obtained with a novel 
canonical algorithm with multiple directed loops, de- 
signed for the calculation of higher-order off-diagonal cor- 
relators. The description of the algorithm is beyond the 
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FIG. 1: (Color online) Superfluid density (upper panel) and 
compressibility (lower panel) for hardcore bosons in a QP 
superlattice; here the system size is L = 1076 and the incom- 
mensuration parameter is a = 830/1076. The cyan horizon- 
tal line marks the critical value V2 / J = 4 for the localization 
transition. The "zebra" -like features in the superfluid density 
are flnite-size effects. 



II. HARDCORE LIMIT 

In this section we present exact results for the hard- 
core case. The study of hardcore bosons in QP potentials 
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has been recently initiated by Refs. 1171181 Remarkably 
the exact mapping onto spinless fermions, Eq. Q, al- 
lows one to exploit the large body of results produced 
several years ago in the context of tight-binding models 
in a QP potentiaP^. In particular, focusing on local- 
ization phenomena, a f undam ental result is the so-called 
Aubry- Andre conjectur^i^E^ stating that for V2 < 4 J all 
single-particle states are extended, while for V2 > 4J all 
states are localized. The single-particle spectrum is or- 
ganized in bands, which have a fractal support?^ (Cantor 
set) at the critical value V2 — 4J. 
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FIG. 2: (Color online) Dispersion relation for a single particle 
in a QP potential of intensity V2/J = 6. 

According to Eq. Q , the properties of the many-body 
system, and in particular its elementary excitations, are 
easily read out from the single-particle spectrum for pe- 
riodic and anti-periodic boundary conditions. Fig. [T] 
shows the superfluid density and compressibility of the 
hardcore-boson system as a function of the chemical po- 
tential and strength of the QP potential. Both quanti- 
ties exhibit a highly non-monotonic dependence on the 
chemical potential /i at all finite values of ¥2- The com- 
pressibility is clearly finite for a finite single-particle den- 
sity of states at the chemical potential, p^^\^) (^^^•'(/i)) 
for periodic (antiperiodic) boundary conditions depend- 
ing on the number of particles as in Eq. [4] in this case 
an extra particle can be added to the system by an in- 
finitesimal change of the chemical potential. As for the 
superfiuid density ps, in a truly periodic system (namely 
in a commensurate superlattice) this quantity is directly 
proportional to the group velocity of the single-particle 
dispersion relation at the chemical potentiapl^, as it im- 
mediately emerges from the definition of ps as response 
function to an infinitesimal phase twist in the operators, 
hi biexp{iSl), namely = (L/2J) d^E/d6^\s=o {E 
is the total energy of the system). Hence ps is finite 
when the chemical potential is inside a band and zero 
otherwis^^. In the case of a QP potential the direct 
connection between superfluid density and group veloc- 
ity breaks down formally, as quasi-momentum is not a 
good quantum number any more, and the single-particle 
states do not have consequently a well-defined group ve- 
locity. Even though the new quantum number k labeling 



the single-particle states is not strictly speaking the mo- 
mentum, the single-particle spectrum e{k) shows a de- 
pendence on k which is reminiscent of that of a system 
in a periodic potential (see Fig. |2] for an example); in 
particular, close to a gap, the dependence of e(fc) on k 
is extremely weak, and it vanishes at the band edge. If 
kcdge is the quantum number of the edge state, and if 
the quantum numbers are defined on an interval A A;, the 
quantum number of the closest state at lower energy will 
^edgc ~ ^odgc ^ Ak/L. A vanishing derivative of the 
energy spectrum with respect to k implies then that 

e(fcedgc) — ^(^cdgc) 

= ^ ^ (7) 

Ak/L 

for L — !■ cx); this means in turn that the energy differ- 
ence e(A:odgc) — ^(^odgo) vanishes faster than 1/L in the 
thermodynamic limit. Consequently, when the chemi- 
cal potential sits close to the band edge, the infinitesi- 
mal perturbation induced by a change in the boundary 
conditions from periodic to anti-periodic is going to mix 
states with an energy difference decreasing faster than 
l/L, and hence is it going to produce an energy change 
^(-) „ which obeys the same scaling. This then 

leads to a vanishing ps as defined in Eq. (|5|. 

On the contrary, for states corresponding to the middle 
of the bands the perturbation induced by a change in the 
boundary conditions will cause an energy shift scaling 
like l/L, leading to a finite Ps, provided that such states 
are affected at all by a local perturbationP^, namely pro- 
vided that they are extended^. Hence in Fig. [l] it is not 
surprising to observe a finite Ps corresponding to a finite 
compressibility for V2/J < 4, and to observe an indenti- 
cally vanishing ps regardless of the chemical potential for 
V2/J > 4, as, according to the Aubry- Andre conjecture, 
the critical point marks a transition from all extended to 
all localized states. This means that the hardcore-boson 
system exhibits an alternation of SF and IBI phases for 
V2/J < 4, while for V2/J > 4 the alternation is between 
IBI phases and BG phases, displaying a finite compress- 
ibility in absence of superfluidity. 



III. SOFTCORE CASE 

When releasing the hardcore constraint for the bosons, 
we lose the possibility of describing the system via free 
fermionic quasi-particles, and a fundamental question 
arises on the fate of the phase alternation described in 
the previous section. For V2 = the phase diagram of 
the one-dimensional Bose-Hubb ard m odel in the grand- 
canonical ensemble is well knowiJSlEl] and it features an 
extended SF region at large J/U ratio, and a succession 
of Mott lobes with integer fillings for lower J/U. Sim- 
ilarly to what is done in absence of a QP potential, we 
study the phase diagram in the (J/U, p,/U) plane, which 
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FIG. 3: (Color online) Density plot for the phase diagram 
of the Bose-Hubbard model in a QP potential with strength 
V2 = U {L = 126). The red line marks the boundary between 
the insulating (small- J/f/) and the superfluid (large- J/17) re- 
gions. The bold black line is the density curve in the atomic 
limit J — > (notice that it is reported on the graph minimum 
J/U = 0.02 for comparison with the numerical data at finite 
J). The dashed blue line shows the same curve for V2 = 
and for a shifted chemical potential ^ ^ fi ~ U/2 (see text). 



is also directly relevant for the intepretation of the behav- 
ior in a trapped system (see Section |v]). For definiteness, 
we choose to study the system with a strong QP poten- 
tial compared to the inter-particle repulsion, namely we 
choose V2 — U in Eq. ([T|. Under this condition the MI 
regions are completely destabilized, given that the Mott 
gap is upper-bounded by U and hence the Mott phase is 
not going to survive the application of an Hamiltonian 
term which is systematically bigger than the gap. Hence 
the only insulating phases surviving in the V2 — U case 
will generically be compressible or incompressible incom- 
mensurate ones, namely insulating phases with incom- 
mensurate filling factors, as shown in Fig. [3] 

We map out the phase diagram for V2 = U via grand- 
canonical quantum Monte Carlo simulations for a single 
lattice size, L = 126 and for a commensurate approx- 
imant a — 97/126 to the incommensuration parameter 
used in the experiment!^. The phase factor of Eq. ^ 
is set to for simplicity. Similarly to the softcore case, 
the characterization of the phases is carried out by inves- 
tigating the compressibility k and the superfluid density 
Ps- The phase diagram emerging from the simulations is 
of extreme complexity, and several observations are in or- 
der, as it will be discussed extensively in the next subsec- 
tion |III A[ It is important to stress that what we present 
in Fig. |4| is the phase diagram for a finite-size system 
and for a single realization of the QP potenti al (nam ely 
a single value of 0); we will see in subsection III A how 



to overcome these limitations. In the phase diagram, in 



FIG. 4: (Color online) Phase diagram of the Bose-Hubbard 
model in a QP potential with strength V2 = U {L = 126). 
Blue dots mark points which are found incompressible. Other 
symbols as in Fig. |3] 



absence of superfluidity, we observe a patchwork of com- 
pressible and incompressible insulating regions. In Fig. |4] 
we indicate with a blue dot a point in parameter space 
which is found to be incompressible. 



A. Finite-size effects 

Before discussing the fundamental features of the 
phase diagram, it is important to stress the following 
technical point. At strictly zero temperature, the num- 
ber of particles is a constant, given that it is a good 
quantum number of the Bose-Hubbard Hamiltonian in 
an external potential; the variation of N upon changing 
the chemical potential fi proceeds in integer steps, and 
hence the compressibility, giving the slope of such curve, 
is simply a succession of (5-peaks. Given a system size L 
and a value such that the ground-state particle number 
is N, the energy gaps between the ground state in the A^- 
particle sector and the ground states in the sectors with 
A^ + 1 and A — 1 particles can either scale to a constant 
for L ^ 00, in which case the system is incompressible 
in the thermodynamic limit, or they can scale to 0, in 
which case the system is compressible in that limit, and 
the apparent incompressibility at flnite L is a flnite-size 
artifact. The Quantum Monte Carlo method we make 
use of is intrinsically a flnite-T approach, and, even for 
a flnite-size system, it can deliver a flnite compressibility 
value if T is larger than the flnite-size gap between differ- 
ent particle number sectors. In the simulation of, e.g., the 
Bose-Hubbard model without QP potential, one can rely 
on the (flnite-temperature) compressibility to approxi- 
mately locate the zero-temperature phase boundary be- 
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tween the compressible SF phase and the incompressible 
MI phase. The error that one can make in that case is 
to ascribe to the SF phase that portion of the MI phase 
whose gap is comparable to the temperature, but this is 
definitely a tolerable error as long as the topology of the 
phase diagram is simple {e.g. one single boundary line 
dividing two extended phases). 

In the case under consideration, we are far from sim- 
pUcity. In fact. Fig. |3] and Fig. |4] are apparently conflict- 
ing, since Fig. |3] seems to suggest an insulating region 
with a smoothly varying n(fj,) over most of the param- 
eter space, namely an insulating region which is mostly 
compressible, whereas Fig.|4]suggests rather the opposite, 
namely that a dominant part of the insulating phase di- 
agram is incompressible. The compressibility data upon 
which Fig. |4] is based are obtained at a temperature 
(T = J/96) which is possibly lower than most finite-size 
gaps in the insulating region, so that the incompressible 
portion of the finite-size phase diagram is overestimated 
with respect to that in the thermodynamic limit. Raising 
T is not recommendable, given that in this way narrow 
incompressible domains in the phase diagram, expected 
in principle on the basis of the hardcore limit (see Section 
[ll| ) , might be washed out if their gap becomes lower than 
the temperature. Hence, given the tight alternation of 
compressible and incompressible regions, the finite-size 
or finite-temperature error that one can make in overes- 
timating the ones over the others is considerable. 



potential, namely the properties of the system in the 
thermodynamic limit become independent of the phase cj). 
This means that increasing the system size towards the 
thermodynamic limit at fixed is equivalent to averaging 
over the phase (/) on a sample of fixed size L, as long as L 
is chosen to be much bigger than any correlation length f 
in the system. Hence, away from the critical line between 
the SF phase and the insulating one(s) (where ^ oo), 
the behavior of the system in the thermodynamic limit 
can be obtained by averaging the finite-size results over 
all possible values of </>. In practice, for a compressible 
region of the phase diagram one expects that, for some 
values of cj), the finite-size gap between different iV-sectors 
becomes smaller than the selected temperature, so that 
a finite compressibility is recovered on average. This nu- 
merical procedure is computationally demanding, and we 
limit it to two selected cuts through the phase diagram, 
one for fixed J/U — 0.04 and one for fixed n/U = 1. The 
results are shown in Fig. |5j they confirm the expectation 
that the fixed-0 phase diagram overestimates the incom- 
pressible regions, but at the same time they confirm the 
presence of sizable incompressible (IBI) regions, in tight 
alternation with the compressible ones (BG). Given the 
finiteness of the grid used in the (J/U,^/U) plane to 
sample the phase diagram, we cannot exclude that an 
even tighter alternation between IBI and BG regions ex- 
ist among the dots, possibly revealed by a finer mesh. 
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FIG. 5: (Color online) Compressibility for the Bose-Hubbard 
model in a QP potential with fixed spatial phase (j> — 0, com- 
pared with the same quantity averaged over tj) fluctuations. 



One possible strategy to circumvent this intrinsic co- 
nundrum of finite-size simulations is based on the follow- 
ing observation. If a truly irrational a is taken in Eq. 
increasing the system size leads to a self-averaging QP 



B. Emergence of the incompressible regions 

The tight alternation of compressible and incompress- 
ible insulating regions upon varying the chemical poten- 
tial reveals that there are several "magic" incommen- 
surate fillings at which a quasi-particle band is filled, 
similarly to what happens to the exact fermionic quasi- 
particles in the hardcore limit. In the compressible BG 
regions, on the other hand, varying the chemical poten- 
tial leads to variations in the fillings, and all the newly 
added particles are effectively localized by the QP poten- 
tial and by the interaction with the density background 
of the other particles. Fig. |3] shows the filling of the 
lattice in the region of phase diagram we sampled nu- 
merically, along with the atomic limit both for the case 
V2 = U and for the case V2 = 0. It is to be noticed that 
the QP potential of Eq. ([2| is symmetric around zero en- 
ergy, and in particular the local chemical potential in its 
minima takes the value /i — U /2; hence a correct com- 
parison with the case V2 = requires to consider a Bose- 
Hubbard Hamiltonian with a shifted chemical potential 
II fi ~ U/2. The presence of the strong QP potential 
alters significantly the density curve in the atomic limit, 
washing out the steps associated with Mott lobes in the 
case V2 — 0, given that the system has lost translational 
invariance and the average density can take all possible 
values. Therefore a seemingly continuous curve emerges, 
with integer densities appearing only at semi-integer val- 
ues of /i/t/; this behavior refiects the broad distribution 
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of local chemical potentials in the sites of the QP poten- 
tial, gradually occupied upon increasing until the n-th 
"shell" is completed when an integer density n is reached. 
The density of atomic states in a QP potential is peaked 
around the energies ±¥2/2 (see Fig.|7]below), as expected 
given that these are the values where the potential has 
minimum derivative, and consequently the density curve 
has maximum slope around fj, — (m + 1/2)U (with m 
integer). Hence the atomic limit picture suggests a fi- 
nite compressibility for all chemical potentials, and intu- 
itively one would expect the compressibility to stay finite 
a fortiori in presence of quantum fluctuations, turned on 
through the hopping term. This would be generally true 
in the case of a truly random potential, on average over 
the possible realizations of the potential or in the ther- 
modynamic limit. Yet, in the case of a pseudo-random 
potential, one of the main quantum effects is to break 
the continuous classical many-body spectrum into bands. 
This mechanism is discussed in the next section. 



IV. TRUE DISORDER VS. PSEUDO-DISORDER 

In this section we describe the formation of bands in 
the quantum many-body spectrum of the Bosc- Hubbard 
model in an incommensurate potential making use of a 
generalization of the atomic limit which includes the first 
quantum correction in a non-perturbative way. The cen- 
tral idea is to subdivide the QP potential in uncorrelated 
wells which are assumed to be independent from each 
other in the case V2/J ^ 1, and to diagonalize the A^- 
boson problem in each of the well. We call this approach 
the random atomic limit. This approach not only repro- 
duces the main features of the n(fi) curve and of the asso- 
ciated many-body density of states at low fillings, but it 
also nicely elucidates the fundamental difference between 
a pseudo-random potential and a truly random one, and 
the fundamental mechanism leading to BG physics. Fi- 
nally it suggests a strategy to emulate the behavior of a 
truly random potential making use of a combination of 
two superimposed QP potentials. 



A. Random atomic limit 

We propose here an approximate treatment of the 
strongly interacting Bosc-Hubbard model in an intense 
QP potential, U — V2 ^ J, m the case of low fillings. 
The idea is that the QP potential subdivides the lattice 
into potential wells between which tunneling is negligi- 
ble when not assisted by the interaction (namely for low 
fillings) and hence inter-well coherence can be neglected 
within a good approximation. Hence the Bose-Hubbard 
model in the grand canonical ensemble can be treated as 
follows: 1) isolate each potential well, cutting the chain 
at the location of the maxima delimiting the well itself, 
and diagonalize the Bose-Hubbard Hamiltonian in this 
potential for TV = 1,2,... particles, finding the ground 
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FIG. 6: (Color online) Schematic depiction ol the random- 
atomic-limit approach for a QP potential (upper panel) and 
for a truly random one (lower panel). This approach subdi- 
vides the system into potential wells, here marked in red, and 
it diagonalizes the problem of 1, 2, and 3 particles confined in 
each well, extracting the corresponding ground-state energies 
-El, E2, and E3. Such energies are then used to reconstruct 
the density of states of the system for low filling. 



state energies Ei, E2, 2) due to the QP potential 
and to the interaction U, adding a second particle to the 
well costs more energy than adding the first, and so for 
the third, namely Ei < E2 ~ Ei < E3 - E2 ~ Ei < ... . 
Hence one can build an effective single-particle density of 
states (DOS) p^^\E) by accumulating the energy costs 
for the addition of a single particle to each well ei = Ei, 
€2 — E2 — El, 63 = E3 — E2 — El, etc. We moreover 
normalize the DOS to unity, namely; 



(8) 



1=1,2,3,. 



The filling for a given fi is simply found by integrating 
p'^^\E) up to p,, similar to the case of non- interacting 
fermions. 

For J = 0, and in the case U = V2 and n < 1, any new 
particle added to the ground state of the system occu- 
pies the least energetic site that remains empty. Hence 
the DOS p'^^^E) in that case simply reflects the distri- 
bution of the on-site energies in the potential, which is 
a continuous one in the interval [—U/2, U/2]. The quan- 
tum corrections to it, obtained for up to TV = 3 parti- 
cles per well, are shown in Fig. [7] the calculation was 
made for a large sample of the QP potential considered 
so far {J/U = 0.04), such that self-averaging of the sys- 
tem properties is guaranteed. We observe the striking 
feature of the opening of gaps in the DOS, defining effec- 
tive bands which are populated by the gradually added 
particles. It is important to stress that the states asso- 
ciated to these bands are localized, and that the nature 
of these bands is determined by the many-body physics 
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FIG. 7: (Color online) Random-atomic-limit result for the 
effective single-particle density of states of a Bose-Hubbard 
model in a QP potential with V2 = U, and over a large system 
size (L = 10''). The quantum case of J/U = 0.04 is compared 
with the classical limit J = 0. 
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FIG. 8: (Color online) Density curve for the Bose-Hubbard 
model in a QP potential with a — 97/126 and V2 ~ U , on 
a system size L — 126. In the inset: random-atomic-limit 
calculation on a big system size (L = 10^). 



inside each well. In particular the first band is associated 
with states in singly occupied wells, and the finite energy 
cost to add an extra particle to any of the well in the sys- 
tem causes the appearence of a gap to the next band. The 
potential-energy cost of adding an extra particle is com- 
pletely compensated for by the gain in chemical potential, 
hence the only remaning energy cost which prevents par- 
ticle addition is purely quantum in nature, and it comes 
from the loss of kinetic energy in the first particle when 
adding a second one to the well. The same reasoning ap- 
plies when adding a third particle, which gives rise to a 
further band, etc. In periodic systems, the emergence of 
Bloch bands is due to the strong scattering of traveling 
matter waves with wavevector at the edges of the Bril- 
louin zone; here we have on the contrary the appearence 
of bands associated with tightly localized particles, and 
band gaps emerging due to strong correlation among the 
particles. 

The opening of gaps in the effective single-particle DOS 
leads to plateaus in the n(fi) curves, and hence to incom- 
pressible phases at special, incommensurate fillings. The 
n(fj.) curve obtained in the random atomic limit is com- 
pared with the QMC results at J/U = 0.04 and for low 
fillings in Fig. |8] We observe that the random-atomic 
limit calculation reproduces very well all the fine details 
of the n(^) curve from QMC: 1) it perfectly captures the 
large plateau which corresponds in Fig. [7] to the gap be- 
tween the 1-particle band and the 2-particle band, occur- 
ring at a filling n = 0.229...; 2) moreover it also correctly 
captures the finite-size plateaus due to the discrete jumps 
in the particle number in the n(ii) curve. Here the sig- 
nificance of finite-size issues in the determination of the 
phase diagram of the system is particularly evident: in 
fact for the L — 126 chain the finite-size gaps are com- 
pletely masking the second plateau in the n(fj,) curve at 
density n = 0.457..., which on the contrary is clearly ex- 



hibited in the random atomic limit calculation on a large 
system size (inset of Fig. |8]). 




E/U 

FIG. 9: (Color online) Density of states from the random 
atomic limit for a truly random potential (upper panel) and 
for a 2-color QP potential (lower panel). The random po- 
tential takes values over the interval [—(7/2,(7/2], while the 
2-color QP potential has V2 = U. For both systems L = 10^. 
Results from two different criteria for the identification of the 
potential wells (see text) are compared with the classical limit 
(J = 0). 
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B. Prom pseudo-disorder to true disorder: the 
two-color quasi-periodic potential 

The discussion of the previous section allows to clearly 
identify the origin of the incompressible behavior in the 
insulating regime of the QP potential. At some specific 
fillings the addition of an extra particle requires to locally 
squeeze the particles already present in the QP potential 
wells, and this squeezing costs a finite amount of kinetic 
energy which translates into a gap over the ground state. 
Actually the energy cost for the addition of a single par- 
ticle can be arbitrarily lowered if the potential wells are 
made arbitrarily large, given that in this way the kinetic 
energy cost for the particles already present can be min- 
imal. In a well of characteristic size I a particle has a 
confinement energy ~ /2mP; the addition of a second 
particle, which is ideally impenetrable to the first one, 
roughly lowers the effective space available to the latter 
to 1/2 and hence it increases its confinement energy to 
^ Afi^ /2'mP. Hence the energy increase scales to as 
when increasing the well size. 

The QP potential has a characteristic length scale 
given by the quasi-period of the QP potential, namely 
I — 211 1 (fcL — fcqp) = (1 — Q^)~^, coming from the beating 
between the incommensurate cosine potential at wavevec- 
tor fcqp = 27ra and the underlying lattice (formally at a 
wavevector fcL = 27r). Potential wells cannot exceed this 
length scale, and therefore the confinement energy scale 
is bounded from below. This means that gaps to par- 
ticle addition at particular incommensurate fillings are 
unavoidable in the system. Moreover, working with a 
fixed number of particles at the same special filling, a 
gap opens to particle-hole excitations: intra-well excita- 
tions are clearly gapped due to the finite size of the well, 
and particle-hole excitations which cause the transfer of 
a particle from a well to another are also gapped, due 
to the chemical potential mismatch between particle and 
hole. 

One fundamental feature of a truly random potential, 
on the other hand, is the absence of an upper bound for 
the extent of potential wells, and hence the possibility of 
always adding a particle to the system at an infinitesimal 
energy cost, or, at fixed filling, to introduce a particle- 
hole excitation at arbitrarily low energy. This feature, 
based on the statistics of rare regions (large w ells), is 
at the heart of the gaplessness of the BG phasd^^^^l^. 
To corroborate this statement quantitatively, we apply 
the random-atomic-limit approach described in the pre- 
vious subsection to a truly random potential with a ran- 
dom on-site energy evenly distributed over the interval 
[— J7/2, U /2\. We use two different criteria to identify the 
potential wells: (criterion 1) a, potential well is delim- 
ited by two successive local maxima with positive height; 
(criterion 2) a potential barrier is identified with a local 
maximum whose height over the previous local minimum 
is larger than C//2, and we hence identify a well with a 
region between two such barriers (see Fig. [6]). Criterion 1 
is less restrictive than criterion 2, and it surely underes- 



timates the size of the potential wells. Yet, in the case of 
the QP potential of Eq. ^ both criteria would lead to the 
same identification of the potential wells as the one used 
in the previous subsection. As shown in Fig. [9] regardless 
the criterion the resulting effective single-particle DOS 
p*^^^ associated with many-body intra-well states is con- 
tinuous for a truly random potential, and all gaps are 
removed by the absence of an upper bound in the well 
length I. 

From the point of view of the cold-atom experiments, 
realizing optically an ideal white-noise potential requires 
to superimpose a large number of standing-wave com- 
ponents with different wavevectors and with the same 
weight; this translates into the need of a large amount of 
lasers at different frequencies and with the same inten- 
sity, which is extremely demanding. Hence an interest- 
ing question, both at the fundamental level and at the 
practical one, is the following: can we mimic de facto the 
physics of a system of bosons in a truly random potential 
by just using a finite number of superimposed standing 
waves, or, in more suggestive terms, how many colors do 
we need for a pseudo-random potential to call it disorder? 

As discussed above, our goal is to realize a pseudo- 
random potential with the minimal amount of Fourier 
components giving rise to a continuous many-body spec- 
trum, namely not exhibiting incompressible IBI phases. 
The guiding principle is that of realizing a potental whose 
spatial features are richer than those of a simple QP 
potential, such that larger wells are realized giving rise 
to intra-well excitations at lower energy, and such that 
states with different number of particles in different wells 
become quasi-degenerate, giving rise to low-energy inter- 
well particle-hole excitations. 

In this perspective we explore the most straightfor- 
ward generalization to the QP potential so far considered, 
namely we add a second incommensurate component to 
the QP potential of Eq. ^ , having the same intensity as 
the first one, thus realizing a 2-color QP potential 

g^'^^i; a, cfi^a' , </>') = cos^(7ra « + 0) +cos^(7ra' i + (t)') - I 

(9) 

In what follows we take a' = 71/126 = 0.56349... 
Fig.[9]shows the random-atomic-limit calculation for such 
a potential on a large sample size (which makes the par- 
ticular values of the phases (p and (p' irrelevant). For both 
criteria of well identification we obtain a gapless single- 
particle DOS, namely no incompressible phase at low fill- 
ing, opposite to the case of a simple QP potential. This 
encourages us to explore the behavior of the system over 
a larger filling interval, making use of QMC simulations. 
Fig. [To] shows the compressibility, superfluid density and 
total density for the Bose-Hubbard model in a 2-color 
QP potential for a large n/U range; the results shown 
are averaged over independent fiuctuations of the 0, (p' 
phase, hence exploiting the whole statistics associated 
with the potential of Eq. ( [9|, and removing finite-size 
effects, as discussed in Sec. |III| Noticeably, the density 
curve at low fillings (n < 0.5) agrees very well with that 
predicted by the random-atomic-limit calculation. Fur- 
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FIG. 10: Compressibility and superfluid density (upper panel) 
and total density (lower panel) for the Bose-Hubbard model 
in a two-color QP potential, Eq. on a L — 126 lattice 
with a = 97/126, a' = 71/126, V2 = U, and J/U = 0.04. 
The QMC results are averaged over independent fluctuations 
of the spatial phases (j), <j>' . 



thermore, over the larger grid of chemical potentials and 
fillings which is accessible to QMC (we went up to n w 3) 
we find a finite compressibility everywhere, generally co- 
existing with insulating behavior (only around fj,/U ~ 2 
a tiny superfluid density shows up). Hence this particu- 
lar pseudo-disordered potential features a continuous BG 
phase without intermediate IBI regions. This means that 
the Bose-Hubbard model in a 2-color QP potential seems 
to capture some of the fundamental physical features of 
the classic dirty-boson problem^. It remains to be seen 
whether the transition from SF to BG in such a potential 
belongs to the same universality class as that of the sys- 
tem in a truly random potential, and whether the absence 
of incompressible regions in the phase diagram applies to 
the whole parameter space: we leave these questions open 
for future investigations. 



V. TRAPPED SYSTEM AND LOCAL-DENSITY 
APPROXIMATION 

In the following sections we discuss the behavior of a 
system of bosons in an incommensurate optical superlat- 
tice in presence of a parabolic trapping potential, to make 
contact with the experiinental setup of current optical 
lattice experiments^^^. Hence we consider the Hamilto- 
nian 



where Vt is the strength of the trapping potential and 
is the center of the trap. 

In absence of an incommensurate potential, the be- 
havior of the Bose-Hubbard model in a trap is generally 
connected to that of the bulk system via a local-density 
approximation (LDA)P2l2Sinil. Here we discuss how that 
approximation holds only partially in presence of an in- 
commensurate potential, and how the behavior of the 
trapped system has to be regarded in its own respect. 
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FIG. 11: (Color online) Local quantities for a Bose-Hubbard 
system in a trap with parameters U = 20J, Vt — 0.0014 U, 
and Htmp — 2.5 U (from top to bottom panel : density, on- 
site particle-number fluctuations, local coherent fraction, and 
local compressibility), compared with the same quantities for 
a bulk system with a chemical potential matching the local 
one in the trap. 



A. Bose-Hubbard model without QP potential 

In presence of a smoothly varying potential, namely 
unde r the assumption that Vt <^ J,U, recent 
works'^ have shown that a local-density approxima- 
tion fully applies to the standard Bose-Hubbard model 
(Eq. ^ with V2 = 0), namely the local properties in 
the trap, as e.g. the density, can be directly related to 
those of the homogeneous system in the grand-canonical 
ensemble by matching the local chemical potential in the 
trap, Hi = /itrap — Vt {i — io)"^, to that of the homogeneous 
system. A similar calculation is shown in Fig. |11[ where 
the local density (n,) and the on-site particle number 
fluctuations 



Htrap = 'HQ + Vt{i- iof 



(10) 



(11) 
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are compared for a trapped system and a bulk system. 
A good agreement between the data for the trapped sys- 
tem and those for the bulk system is observed away from 
the crossover regions between a locally SF and a locally 
MI region, while in the crossover regions the data for the 
trapped system have generally less sharp features with 
respect to those for the homogeneous one. This can be 
intuitively understood by noticing that finite-size effects 
are typically most pronounced close to critical bound- 
aries between different phases, where discontinuities and 
singularities are rounded off by the finiteness of the sys- 
tem. In the trapped system the finite-size effects are not 
only given by the overall trapping potential, but most 
significantly by the finiteness of those regions whose lo- 
cal chemical potential corresponds to critical boundaries 
between MI and SF in the homogeneous case. 

Despite the finite-size corrections, the above results 
lead typically to interpret the various regions in the 
trapped system as locally superfluid or locally Mott in- 
sulating depending on the corresponding phase at the 
same chemical potential in the phase diagram of the bulk 
system. This conclusion is particularly confirmed when 
looking at the local condensate fraction 



(12) 



and at the local compressibility, introduced in Ref. 1301 

13 J ■ 



'^loc,! 



L 



^{{n.rij) - {n,){n.j)). (13) 



particle-hole excitations, and it is typically accompanied 
by a lower local phase coherence. These two combined 
observations form the basis to the LDA interpretation of 
local MI and SF phases in trap. This interpretation is 
fundamental to conclude that the MI phase is realized at 
all in the trapped system, with a halo of SF surrounding 
it; indeed only an infinite repulsion U and an infinitely 
steep trap would allow to realize a pure MI in the trap. 
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Notice that in the homogeneous bulk system with 
global condensate fraction n{k = 0) and global compress- 
ibility K, we have that niacAk = 0) = n{k = 0) / L and 
Kioc.z = n/L. From Fig. [ft] we see that the behavior 
of both n\oc.i{k — 0) and n\oc.i in the trapped system 
follows qualitatively that of the global quantities in the 
bulk system. In particular it is important to notice that 
both nioc.i(fc = 0) and Kioc.i are actually non-local quan- 
tities, and their final values depend not only on the lo- 
cal chemical potential at point i but also on the prop- 
erties at other points in the trap, and in particular on 
the total number of particles in the trap. Thus a close 
agreement with the bulk data is generally not expected, 
in particular whenever the correlation functions {b\bj) 
and {{niTLj) — {ni){nj)) become long-ranged in the bulk 
system, namely in the SF phase and at the SF-to-MI 
boundary, respectively. Nonetheless nioc,i(fc = 0) and 
Kioc,i capture the short-range behavior of those correla- 
tors in the trapped system around the point i; therefore 
their local enhancement expresses the fact that phase co- 
herence between the points in those regions is strong, 
which is typical of a locally superfluid islands, and that 
particle number fluctuations in those regions are strongly 
correlated, which is a signature of the existence of low- 
energy and long- wavelength density modes. Conversely 
a vanishing local compressibility manifests a local gap to 



FIG. 12: (Color online) Local density, particle-number fluctu- 
ations and coherent fraction (from top to bottom) for a Bose- 
Hubbard system in a trap plus QP potential, with parameters 



U = 30J, Vt = 0.0014 U, 



1.8 U, and V2 = U. The 



above quantities are compared to those of the bulk system 
with matching local potential (Eq. (14|). 



B. Bose-Hubbard model with QP potential 

When introducing the incommensurate cosine poten- 
tial as in Eq. ([T|) , the picture becomes significantly more 
complex. Wc will henceforth limit our discussion to the 



case V2 — U already discussed in Sections III and IV 



namely the case in which the MI phase is completely re- 
moved from the phase diagram and substituted by a va- 
riety of incommensurate gapped and gapless insulating 
phases. 

In presence of an incommensurate cosine potential, 
even the bulk system becomes inhomogeneous, and hence 
the LDA has to be rephrased as the approximation re- 
lating the properties of a point in the trapped system to 
those of a point in a bulk system experiencing the same 
local chemical potential, resulting from the global chem- 
ical potential plus the cosine potential; in mathematical 
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FIG. 13: (Color online) Upper panel: local density of the 
trapped system compared with local density and local com- 
pressibility of a bulk system with chemical potential /Xbuik = 
1.8 U. The region where trapped and bulk density profiles 
match is emphasized by the shaded area (notice that the 
chemical potentials /itrap and /^buik match perfectly only at 
the trap center, but the variation of /xtrap due to the trapping 
potential is slow around the center). The regions exhibiting a 
finite local compressibility in the bulk system are seen to fall 
out of the region of matching densities, and consequently the 
trapped system is locally incompressible (see following panel). 
Lower panel: local compressibility of the trapped system com- 
pared with the global compressibility of the bulk system. All 
other parameters as in Fig. |12[ 



terms, we have to identify a site itrap of the trapped sys- 
tem with a site ibuik of the bulk system such that 

Woe = A^trap + Vt (itrap - iof + V2 g{itra.p] O., (j)) 

= A^buik + V2 g{ihuik; a, cj)) (14) 

where ^tiap ^^nd /ibuik ^^re the global chemical potentials 
for the trapped and the bulk system, respectively. Fig.[T2| 
compares the results for the bulk and the trapped system 
following the prescription of Eq. ( 14 1 in the case of a 
strongly interacting system in a dc 
potential, U/J = V2/J = 30. 



; eep incommensurate 
With these parameter 
values the bulk system lies in the insulating region for 
a broad range of chemical potentials. In what follows 
we consider /.ttrap/f^ — 1-8 so that the local chemical 
potential in the trapped system, which is upper-bounded 
by the above value, experiences a superfluid region of 
the bulk phase diagram (Fig. |4]) only around the value 
Woc/C^ = 1-5. 

Despite the high inhomogeneity of both the bulk and 



the trapped system, we generally observe a good quan- 
titative agreement between their respective data for the 
local average density, local density fluctuations and local 
coherent fraction. For the interval of chemical poten- 
tials spanned by /iioc only two nearby critical points are 
present in the bulk phase diagram around fi\oc/U = 1.5, 
so that, away from those two points, the system is deep 
in an insulating region with a short correlation length, 
and finite-size effects are therefore not severe. The prob- 
lematic aspect of LDA in presence of the QP potential is 
nonetheless the subsequent interpretation of the trapped 
system as locally exhibiting a phase of the bulk system. 
In particular, what is the meaning of a local Bose glass! 

In this respect, the examination of the local compress- 
ibility compared with the global one is illuminating. As 
shown in Fig. |13| the local compressibility of the trapped 
system is in strong disagreement with the global one, and 
it is generally vanishing when the global one is not. This 
is not at all surprising in a strongly inhomogeneous sys- 
tem, in which the global compressibility is due to the 
local response of disconnected regions to the variation of 
the chemical potential, and it is related to the existence of 
localized low-energy particle-hole excitations. Although 
a region in the trap might be at a chemical potential 
which would correspond to a BG phase in the bulk phase 
diagram, it is very likely that such region does not cor- 
respond to the region of the bulk system that is hosting 
the locally quasi-gapless excitation and which is hence ex- 
hibiting a finite local compressibility. This is indeed the 
case for what depicted in Fig. [13] where we show the local 
compressibility of the bulk system for a chemical poten- 
tial corresponding in the trap to a locally incompressible 
region: as observed there, the bulk system has a finite 
compressibility at that chemical potential, coming from 
the local compressibility of regions which do not corre- 
spond to the ones reproduced in the trap (shaded areas) . 
Hence it is evident how the finite-size effects induced by 
the trap are drastically altering the local behavior of the 
system with respect to that of the bulk one. 



C. Structure of the particIe-hoIe excitations 

The above results require to discard the LDA picture 
of the trapped system locally mimicking the behavior of 
the bulk system when it comes to the BG phase. At 
the same time, the very picture of a local realisation of 
collective quantum phases in a trapped system might be 
a limiting point of view. As pointed out in Section IV 



the existence of gapless particle-hole excitations originat- 
ing from rare regions in the system is at the core of the 
BG phase: e.g. accidentally degenerate regions sepa- 
rated by a potential barrier, for which tunneling creates 
a pair of quasi-degenerate symmetric/ antisymmetric su- 
perposition states similarly to the double-well problem; 
or rare regions with a locally uniform or periodic poten- 
tial, which host local low-energy excitations. 

As we have seen, the inhomogeneity of the QP poten- 
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FIG. 14: (Color online) Central panels: quasi-particle 

(IV'quali-p(«)n and quasi-hole (-lV'quii-h(OP) densities for 
the Bose-Hubbard model in a trap at various particle numbers 
A'', and with parameters U = 20J, Vt = 0.0014 U. A rescaled 
picture of the parabolic trap is reported in each panel for 
reference. Upper and lower panels: total density profile asso- 
ciated with the minimum (A'' — 51) and maximum (A^ = 62) 
number of particles shown here. 




FIG. 15: Particle-hole gap Ai5p_h for the Bose-Hubbard 
model in a trap. Parameters as in Fig. |14| 



tial is what causes a trapped system to fail in reproduc- 
ing locally the global properties of a bulk system at the 
same chemical potential. At the same time it is conceiv- 
able that, in presence of the trap, accidental degeneracies 
might occur between local quasi-particle states living at 
different trapping potentials, namely accidental degen- 
eracies which are not observed in the bulk system. To di- 
rectly investigate this possibility, we have run a canonical 
Monte Carlo simulation in which we add particles pro- 



gressively to the system. We then define effective quasi- 
particle and quasi-hole states ipl^Li-p ^^'^ ''P'^Li-h such 
that their associated density distribution corresponds to 
the density variation of the system by adding/removing 
a particle, respectively: 



^^^^ H) 

V'iuali_h(«) 



(AT-l-l) (AT) 

Ti: — rii 



(N) (N-1) 

n) — n] 



(15) 
(16) 



Looking at the spatial structure of the quasi-particle and 
quasi-hole state corresponding to the same particle num- 
ber TV, we generally obtain information about the struc- 
ture of the particle-hole excitations of the system at par- 
ticle number N . 

In fact, starting from iV — 1 particles, the energy cost 
of adding an extra particle is given by the kinetic en- 
ergy plus the interaction of the extra particle with the 
background of the — 1 particles already present in the 
system, and analogously for adding an extra particle to 
go from A^ to A^ -I- 1. If the density variation associ- 
ated with growing the particle number from A'^ — 1 to 
A'^ happens in a different spatial region with respect to 
the variation caused by growing the number from A^ to 
A^ + 1, then the quasi-hole state V'q^yisi-h quasi- 



particle state ip 



(N) 

quasi— h 



are non-overlapping. In pres- 



ence of short-range interactions as in the Bose-Hubbard 
model, this implies that the addition of the A^-th par- 
ticle does not influence the background with which the 
(A^ + l)-th particle interact, which implies, in more co- 
incise terms, that the quasi-particle and the quasi-hole 
are not interacting. Thus, working at particle number 
N, the creation of a quasi-hole in state V'q^isi-h ^'^^ ^ 

quasi-particle in state V^q^asi-p represents the true lowest- 
energy particle-hole excitation of the system, with energy 
^En = A^p_h = S^li + Ef_^ - 2eP. 

If, on the contrary, the quasi-particle and quasi-hole 
states for particle number A'^ do overlap, then in prin- 
ciple the energy Aii'p_ij needs a correction coming form 
the interaction in order to be identified with the energy 
of the elementary particle-hole excitation. Yet, if A_Ep_h 
is found to be small (<C U for definiteness) , then the 
interaction between the A^-th added particle with the 

background of A^ — 1 particles, entering in i?^' — -E^^^ 
difference, has essentially the same energy as that of the 
interaction between the (A^ + l)-th particle and the back- 
ground of A^ particles, entering in — E'j^^ difference; 
this in turn implies that the A^-th particle minimally al- 
ters the interaction between the (A^ + l)-th particle and 
the remaining A^ — 1, namely the quasi-particle and the 
quasi-hole states are again effectively independent, and 
they are moreover essentially degenerate. If not even 
the above condition on Ai?p_h is satisfied, we can still 
generally assume that, in presence of repulsive interac- 
tions, Ai?p_h represents an upper bound to Ai^Tvi^and 
that the spatial structure of the true particle-hole exci- 
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tations is approximately reproduced by that of 



(N) 



" quasi— h 

and V'quii-p- 

In the case of a trapped system without QP poten- 
tial, the quasi-particle and quasi-hole states are shown 
in Fig. [14] for a selected range of particle numbers N, 
and the associated particle-hole gap Ai?p_h is shown 
in Fig. [15] We generally observe that the quasi- 
particle/quasi-hole states live in the locally superfluid 
regions, as expected from the information on the local 
compressibility. For the particular density we have cho- 
sen in Fig. 14 the system exhibits two such regions, at 



local filling n < 1 on the wings and at filling n > 1 in 
the center. We generally observe three different types 
of particle-hole pairs, with AE'p-h <C t/ in all cases: 
1) pairs where quasi-particle and quasi-hole both live 
in the center of the trap — these are the excitations 
which more closely mimic those of a finite-size SF sys- 
tem; 2) pairs where quasi-particle and quasi-hole both 
live in the wings of the trap — these states correspond 
to the symmetric and antisymmetric superpositions of 
the left- and right-localized states in the wings, and the 
particle-hole gap Ai?p_h is minimal in this case, because 
the large barrier provided by the intermediate MI region 
leads to quasi-degeneracy (see for instance the case of 
N = 55,56,61,62); 3) particle-hole pairs with a particle 
living in the wings and the hole in the center or viceversa. 
The excitations 2) and 3) are specific of the trapped sys- 
tem; in particular it is noteworthy that quasi-degeneracy 
exists between states related to different trap regions, 
hence the picture of local phases described in the previ- 
ous paragraph does not extend to excitations. 

In presence of a strong QP potential with V2 — U, the 
structure of the excitations changes drastically, although 
some of the features observed before persist, namely 
the cross-talk between different trap regions. Fig. [161 
shows the spatial structure of quasi-particle/quasi-hole 
states, and Fig. [iT] the particle-hole gap. It is evident 
that the strong QP potential leads to tight localiza- 
tion of the added particles/holes in the system, typically 
around a few lattice sites; delocalization of the quasi- 
particle/quasi-hole state onto two disconnected regions 
is due to the symmetric-antisymmetric combination of 
quasi-degenerate localized states (as before, these combi- 
nations give rise to the lowest particle-hole gaps, as e.g. 
for N = 69, 72, 76, etc.). Furthermore, quasi-particle and 
quasi-hole states are very often decoupled spatially (see 
for instance the cases N = 70, 71, 73, 74, 75, 77, 78). This 
means that the lowest-energy particle-hole excitations are 
not at all confined to equipotential regions, but they can 
be associated with accidental quasi-degeneracies of quasi- 
particle and quasi-hole states between regions at quite 
different trapping potentials. Hence the notion of "local 
( quasi- )gaplessness" is too restrictive, and the picture of 
a local BG, which we have seen to break down in the pre- 
vious discussion, can be superseded by the possibility of 
having a global trapped BG, in which low-energy particle- 
hole excitations, albeit corresponding to the transfer of 
a quasi-particle between two localized states, can also 
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FIG. 16: Central panels: quasi-particle (| V'q^Lsi-p(*)l'^) ^^"^ 
quasi-hole (— I V'q!L'si-h(*)l^) densities for the Bose-Hubbard 
model in a trap plus QP potential at various particle num- 
bers A^, with V2 ~ U. Other parameters as in Fig. |14[ A 
rescaled picture of the parabolic trap plus the QP potential 
is reported in each panel for reference, together with the trap 
axis. Upper and lower panels: total density profile associated 
with the minimum {N = 69) and maximum (TV = 80) number 
of particles shown here. 




FIG. 17: Particle-hole gap AEp-h for the Bose-Hubbard 
model in a trap plus QP potential. Parameters as in Fig. |16| 



occur between distant points in the trap. 



D. Bose glass in the trap: 1-color vs. 2-color QP 
potential 

In Section |V| we have seen that for a single realization 
of the QP potential the concept local BG (in the sense 
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FIG. 18: (Color online) Upper panel: phase-averaged local 
compressibility for the Bose-Hubbard model in a QP po- 
tential, with parameters U = 25J, /itrap = 1.8 U, and 
Vt = 0.0014 U. Lower panel: the phase-averaged local 
compressibility of the trapped system is compared to the 
phase-averaged global compressibility of the bulk system (see 
also Fig. [5| by matching the average chemical potentials: 
here /left = fJ-tis.p — Vt{i — io)^ for the trapped system and 
Mcff ~ Mbuik for the bulk system. 



of the LDA) breaks down due to significant finite-size ef- 
fects. This feature is actually generic of inhomogeneous 
systems, and we expect it to hold also for a single realiza- 
tion of a truly random potential. The concept of local BG 
can be nonetheless recovered on average over the statis- 
tics of the (pseudo)-disordered potential. In the specific 
case of a QP potential, one can randomize its spatial 



phase as in Sec. IV Interestingly, from the point of view 
of optical-lattice experiments, averaging over fluctuations 
of the spatial phase of the QP potential is inevitable when 
averaging the results over different experimental runs. In 
fact the spatial phase can be fixed over the duration of a 
single run, but it is typically changing from run to rurP^ 

Randomizing the spatial phase of the QP potential 
implies that the potential profile can randomly "slide" 
along the trap. If a region of the trap experiences a lo- 
cal chemical potential which corresponds to that of a BG 
phase in the bulk phase diagram, there is in principle 
a finite probability that, by random sliding of the QP 
potential, a region hosting a gapless particle-hole excita- 
tion in the bulk system will appear in the trap, so that 
on average that portion of the trap acquires a finite lo- 
cal compressibility. Fig. [19] compares the data, averaged 
over (/)-fluctuations, of the rescaled global compressibility 
k/L for the bulk system (already shown in Fig. [5]) with 
those of the local compressibility in the trap, plotted as 
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FIG. 19: Left panel: local compressibility of the Bose- 
Hubbard model in a trap plus a 2-color QP potential, averaged 
over fluctuations of the spatial phases of the two colors. Here 
Vt = 0.0014 U, other parameters as in Fig. |10[ Right panel: 
momentum distribution for the same system. 



a function of the effective chemical potential induced by 
the trap, ficS — Mtrap — Vt{i — jq)^. This time we ob- 
serve that, upon (/)-averaging, most compressible phases 
in the bulk are generally mirrored by locally compressible 
phases in the trap, and even that there exist locally com- 
pressible phases in the trap which correspond to incom- 
pressible phases in the bulk. This suggests that the finite 
compressibility in this case is associated with low-energy 
excitations which are specific of the trapped system, as 
discussed in the previous section. 

The above results finally suggest the following conclu- 
sion: a truly random potential or a pseudo-random po- 
tential, giving rise on average to a BG phase for a contin- 
uous range of chemical potentials in the phase diagram 
of the bulk system, should also be able to give rise on av- 
erage to a locally compressible phase on all the points of 
the trapped system. To verify this statement in a partic- 
ular case, we look at the 2-color QP potential of Sec. |IV| 
in a trap. This potential is seen to satisfy the above con- 
dition on the compressibility over an extended range of 
chemical potentials (see Fig. 



10 1. Fig. 19 shows the lo- 



cal compressibility of this system, averaged over spatial- 
phase fiuctuations. In contrast with what observed in 
Fig. [12] and Fig. [18) we have here a finite compressibility 
everywhere in the trap, on average over (/)-fluctuations. 
At the same time the system exhibits a very low co- 
herence which allows to rule out quasi-condensatiorP^l. 
Hence we conclude that this system perfectly realizes a 
BG all over the trap. 



VI. RELEVANT EXPERIMENTAL 
OBSERVABLES FOR THE TRAPPED SYSTEM 

In this section we show results for two relevant ob- 
servables that are currently accessible to time-of-flight 
measure ments in optical lattice experiments, namely the 
first-ordei *^°*^^l and second-ordei'^S'^MiiiSj coherence of 
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the bosons. The first-order coherence contains the in- 
formation on the momentum distribution, Eq. (|6|. The 
second order coherence contains instead the information 
on the correlations between the momentum distribution 
fluctuations (also called noise correlations): 

G{k,k') = {n{k)n{k')) - {n{k)){n{k')). (17) 

What is costumarily measured in experiment^^HHl is ac- 
tually the average correlations between the populations 
at two momenta differing by a value q, normalized to the 
factorized correlator: 



feSFBZ 



[n{k)n{k + q)) 



E 



feGFBZ 



{nik)){n{k + q)) 



(18) 



In order to evaluate G{k, k') and I{q) we need in gen- 
eral to be able to evaluate 4-point off-diagonal correla- 
tors of the type {b\bjb\b„i) within QMC, and to succes- 
sively Fourier transform them. This is possible within a 
double- directed-loop canonical algorithm that we explic- 
itly developed for this purpose, in which the evaluation 
of the above cited 4-point correlator is performed during 
the update whenever the 4 ends of the two worms find 
themselves at the same time slice on the sites i,j,l,m. 
Details of the algorithm will be reported elsewhere. We 
investigate the evolution of the first- and second-order co- 
herence upon increasing the strength of the incommen- 
surate cosine potential. We fix the ratio U/J to two 
different values U/J = 5, 10, which give rise to two dif- 
ferent states at zero QP potential in a trap of strength 
Vt/J = 0.0014 U. 

For U/J = 5 the trapped system with V2 — exhibits 
a high coherent fraction, and the application of the QP 
potential has the effect of suppressing coherence, driv- 
ing the system through a superfluid-to-incommensurate- 
insulator crossover, as clearly shown in Fig. 20 We 
then compare the evolution of the momentum distribu- 
tion with that of the density profile upon increasing V2 
(also shown in Fig. 20 1, which is not directly measurable 
in the current experiments. We notice that the strongest 
suppression of the coherent peak at A; = corresponds 
roughly to the value of V2 at which the incommensurate 
potential introduces unoccupied sites in the central re- 
gion of the trap, fragmenting the many-body state into 
droplets. 

Along the superfluid-to-insulator crossover the system 
traverses a modulated superfluid phase, in which the co- 
herent fraction remains significant while the density pro- 
file acquires a strong modulation due to the external po- 
tential. This modulation is only marginally revealed in 
the momentum distribution by two satellite peaks at the 
incommensurate beating momenta fcjnc = ±27r(l — a) 
(see also Ref. "STfor a similar signature in commensurate 
superlatticcs) . This effect can be understood as an effec- 
tive reduction of the first Brillouin zone, due to the longer 
(quasi-)period imposed by the incommensurate potential; 
yet its signature is possibly too weak to be observed in 
current optical lattice experiments. It is also to be ob- 
served that such signature becomes unobservable in the 
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FIG. 20: (Color online) Upper panel: momentum distribution 
for the Bose-Hubbard model in a trap (A'' — 100, U = 5J, 
Vt = 0.0014 U) in a QP potential of increasing strength V2 
(here represented in units of J). Lower panel: false-color 
density plot for the same model; the false colors only apply to 
finite densities, while the white regions emphasize the points 
with zero density. 



strongly insulating regime, where the momentum distri- 
bution becomes essentially featureless. 

Nonetheless, as shown experimentally in Refs. I25|45} 
the fluctuations of the momentum distribution, captured 
by Eqs. (17 1,(18 1, have a momentum structure which 



reflects directly the Fourier transform of the density- 
density correlations. Indeed, for a pure Fock state |{?^^}) 
one can easily show that 



N 



N 
L2 



S{q) 



(19) 



liq) = 1 

where S{q) = l/LJ2^^{ 
tor, and N is the total number of particles. In a more 



is the static structure fac- 



general superposition state j^*) : 
the above relation does not hold; 



= E{„,} c({nj) |{nj) 
nonetheless, if the state 
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contains only one or a few dominant Fock components, 
the essential features of their related structure factor will 
be captured by the integrated noise correlations as in 
Eq. (pi. 



Fig. 



21 shows the second-order coherence as 



a function of the intensity of the incommensurate po- 
tential. The central peak at g = shows a non-trivial 
evolution, with a large increase corresponding to the 
suppression of coherence, and the consequent increase 
in the fluctuations at all momenta, which con- 

tribute to I{q = 0). But the most significant feature is 
the appearence of satellite peaks a,t q — hmc and also at 
q = 2A:inc. The height of these peaks appears to saturate 
around the value of V2 (V2/J ~ 20) at which the coher- 
ence peak in the momentum distribution gets drastically 
suppressed, marking the fragmentation of the system into 
droplets. 

Increasing the value of U/J to 10, the system at 
V2 = has a significantly lower coherence than in the 
case U/J = 5. Interestingly the application of the QP 
potential appears to increase the coherence in the sys- 



tem, as shown in Fig. 22 an effect due to the simple fact 



that MI regions are destabilized by the application of the 
QP potential because they obviously do not display the 
right filling to minimize the potential energy. A similar 
effect of destabilization of the MI by an externally ap- 
plied potential has been previously observed both in the 
case of a periodic potentiaP^ and in the case of a ran- 
dom on^^. A further increase in the QP lattice leads 
to a similar crossover towards a deep incommensurate 
insulating state, as seen above for the case U/J = 5. 
The second-order coherence shown in Fig. [23] also reveals 
a non-monotonic behavior of the q — peak, where a 
suppression of momentum fluctuations are observed to- 
gether with the enhancement of the fc = peak. Satellite 
peaks at k^c and 2/cijjc appear for large V2 as in the case 
of a lower U /J ratio, but it is interesting to notice that 



the momentum structure of I{q) is signiflcantly less pro- 
nounced than in the case of U/ J = 5, in particular the 
height of the satellite peaks has decreased with respect to 
that of the central one. To further illustrate this trend we 
have performed an exact calculation of I{q) in the Fock- 
state limit J — > 0, whose results are shown in Fig. [23] 
and where it appears that in this limit the g-dependent 
structure of I{q) is even weaker than in the finite- J case. 
Hence we can conclude that quantum fluctuations gener- 
ally enhance the non-trivial g-dependent part of the I(q) 
signal with respect to the q — peak. 

In conclusion, we observe that the first- and second- 
order coherence provide full information about the di- 
agonal and off-diagonal correlations in the system, and 
they clearly reveal the crossover (or finite-size transition) 
from a superfiuid state to an incommensurate insulating 
state upon increasing the QP potential. The detection 
of the possibly unconventional nature of the insulating 
state, namely the presence (albeit local) of BG regions, 
is nonetheless beyond the scope of these observables, and 
quantities that directly probe the low-energy particle- 
hole excitations are required, ideally the local compress- 
ibility as discussed in Section [V] Future work will present 
an experimental strategy to directly address this issud^. 



VII. VALIDITY OF THE BOSE-HUBBARD 
MODEL IN AN INCOMMENSURATE 
POTENTIAL 

Throughout the previous sections we have worked un- 
der the assumption that an incommensurate optical su- 
perlattice, realized in the experiments^Sl by superimpos- 
ing a primary and a secondary standing wave with incom- 
mensurate wavelength relation, can be fully described by 
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FIG. 23: Upper panel: noise correlations for the Bose- 
Hubbard model in a QP potential with U — lOJ; other pa- 
rameter as in Fig. |22[ Lower panel: noise correlations for the 
classical limit of the above model, J — 0. Here V2 is in units 
of U/W. 



a simple incommensurate cosine potential added to the 
Bose-Hubbard Hamiltonian, as in Eq. This means 
that the only effect of the secondary optical lattice is as- 
sumed to be a shift of the local energy of the Wannier 
functions associated with the lowest Hubbard band of 
the primary lattice. Yet, the frequency inside each po- 
tential well of the superlattice changes in general with re- 
spect to that of the single-color lattice, and, most impor- 
tantly, the relative distances between consecutive wells 
are also shiftecP®. These effects modify in principle the 
shape of the Wannier functions, leading to a modulation 
of the U parameter, and the overlap between the Wan- 
nier functions associated with two adjacent sites, leading 
to a modulation of the J parameter. 

Here we recall the well-known derivation^ of the Bose- 
Hubbard model from the most general second-quantized 
Hamiltonian of a set of interacting bosons in an external 
potential, which reads 



+ g j d^r ijj\r)'ip\r)il^{r)i^{r) (20) 

where V'jV'^ are bosonic field operators, g — inasfi^ /m 
with ni the mass of the bosons and the s-wave scatter- 
ing length, and Vopt{x) is the optical potential applied to 
the atoms. Following the experimental setup of Fallani 
et al^, hereafter we assume that Vopt can be written as 



Kpt(r) = Vii{x) + V^iy,z) 

V±iy,z) — s_L [cos^(fci?/) + COS^(fciz)] ii^j. 

+ ^muj'i [{y - yo)'^ + {z - zof] 

V\\{x) — Si COS^ (kix) Er + S2 C0S^{k2X + ' 



(21) 
(22) 

Er 



V± creates a strong optical lattice potential with 
wavevector ki = 2tt/\i in the y and z direction, defin- 
ing tubes along the x direction, plus an overall parabolic 
confinement with frequency lo_\_. Hereafter we consider 
a strong transverse optica l lattice, s± = 40, as obtained 
in recent experimenta^^'^Il which allows us to neglect the 
inter-tube hopping and to consider single tubes indepen- 
dently. The optical potential along the x direction fea- 
tures two optical lattices with incommensurate wavevec- 
tors fci = 2ir/\i and fc2 = 27r/A2, plus a parabolic con- 
finement with frequency cl;||. For definiteness we take 
Ai = 830 nm, A2 = 1076 nm, and mw2/2 = 0.0012£;r as 
in Ref. ilOi si, S2 and sj_ are dimensionless amplitudes, 
and Ej. = {KkiY /2m. is the recoil energy (associated with 
the primary lattice). 

We then decompose the bosonic field operators onto an 
orthonormal basis of single particle wavefunctions local- 
ized around the minima of the optical potential in each 
tube: 



V'(r-) = ^w;,(a;)M^,(y)l¥,(z) h 



(23) 



where i runs over the minima. Neglecting the overlap be- 
tween localized wavefunctions on non-adjacent minima, 
the Hamiltonian parameters of a general one-dimensional 
Bose-Hubbard Hamiltonian 

7i = ^ (jj.i+ifoifoj+i + h.c.j + Y ni{ni -l) + Vi n. 

(24) 



can be obtained as usual, with parameters 

h"^ d 



Ji.i+i = I dx w*+i(x) 
dx \wi{x)\'^ 



w,{x) (25) 



2ma^ + ^|("). 
^ dy miytj dz m{zt {26) 

= I dx \w,{x)\''V\\ix) (27) 
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FIG. 24: Hopping amplitude J as a function of the lattice po- 
tential depth si for a monochromatic optical lattice (both 
quantities are expressed in units of the recoil energy Er). 
We show the estimate of J obtained through the solution 



exp( 



of the Mathieu equatiorPl, J/Er = (4/^ 
through the calculation of the overlap integral Eq. 
e^'-order expansion of the optical potential (see text 
through the simple Gaussian approximation for the Wannier 
functions, J/Er = (7r^/4 — 1) si exp[— (7r^/4)^/s7]. 
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FIG. 25: (Color online) Site-dependent hopping amplitude 
Ji,i+i and local potential Vi for the generalized Bose-Hubbard 
model, Eq. (24 1, derived from the microscopic Hamiltonian 



for bosons in an incommensurate optical superlattices. Here 
si = 8, S2 = 3 and mc^jj/2 = 0.0012^;^. 



The Wi's are the Wannier functions associated to the 
lowest band of the periodic potential in the y and z direc- 
tions, which can be conveniently approximated by Gaus- 
sians in the limit of a deep optical lattice. Along the 
X direction Bloch's theorem does not apply due to the 
incommensuration between the two optical lattice com- 
ponents, so that Wannier functions cannot be properly 
defined. Hence we take the Wi to be the ground state of 
the potential obtained by a polynomial expansion (up to 
gth Qj.(jgj.^ of V| I (x) around the i-th minimum, and the 



FIG. 26; Upper panel: relative fluctuations of the effective 
hopping amplitude Ji,i+i for bosons in an incommensurate 
superlattices, as a function of the intensities si and S2 of the 
two superimposed standing waves. Lower panel: ratio of the 
hopping amplitude fluctuations over the potential energy of 
the secondary standing wave. 



so-obtained set of localized functions is further Gram- 
Schmidt orthogonalized. In absence of the A2 optical 
lattice, this approach provides hopping amplitudes Ji^i+i 
in good agreement with the estimate coming from the 
solution of the Mathieu equatiorP^l^ and it improves sig- 
nificantly over a simple Gaussian approximation for the 
Wi's (see Fig. 24). In the case of the full incommensurate 



super lattice, this approach shows a significant site depen- 
dence of the hopping amplitudes Ji,i+i, with differences 
between different pairs of sites that can reach one order 
of magnitude for intense secondary lattices. This is due 
to a significant shift in the positions of the minima of the 
incommensurate optical superlattice with respect to the 
monochromatic lattice, and to the exponential sensitiv- 
ity of the hopping amplitude to such shifts. An example 
of the site-dependence of the hopping amplitudes for an 
intense secondary lattice is shown in Fig. 25 together 
with the local energies Vi. In particular, it is evident 
that the modulation of the local chemical potential and 
of the local hopping are strongly correlated, namely the 
amplitude of the hopping is maximum where the local 
energy Vi is also maximum. Hence the modulated hop- 
ping competes in principle with the localization effects 
induced by the QP Vi's, given that the local kinetic en- 
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ergy is minimized in proximity of the potential maxima. 
Yet, looking at the magnitude of the hopping modulation 
with respect to the dominant local potential modulation, 
it is clear that one can still confidently rely on the simple 
Hamiltonian Eq. ([ij to capture the dominant low-energy 
features of the Hamiltonian Eq. (20 1. Fig. 26 shows the 
relative fluctuations bJj Jav associated with the modu- 
lation of Ji.i+i around the average value Jav, and their 
amphtude 6 J compared with the energy scale of the po- 
tential created by the secondary lattice, S2-Er- Although 
the relative fluctuations can be extremely big and far 
exceed 100%, the energy range spanned by those fluctua- 
tions becomes negligible with respect to the energy scale 
of the potential, and hence the modulation of the hop- 
ping is not expected to alter significantly the behavior of 
the system for strong secondary lattices. As for the site 
dependence of the interaction Ui, we find that it reaches 
~ 10% for the strongest secondary lattice considered in 
Fig. |26l so we discard it for simplicity in the following. 





FIG. 28: (Color online) Upper panel: noise correlations for 
bosons in an incommensurate superlattice, same parameters 
as in Fig. [27| (upp er p anel). Lower panel: same quantity 
for the model Eq. ( |24[ ) without accounting for the hopping 
modulation. 




FIG. 27: (Color online) Upper panel: momentum distribu- 
tion of A'' = 100 bosons in an incommensurate superlat- 
tice with si = 8 and variable S2 {s± = 40Er, as = lOOao 
as in *^Rb with ao ~ Bohr radius, and other parameters 



as in Fig. 25 I , ca lcul ated for the generalized Bose-Hubbard 
Hamiltonian Eq. ( 24 1 . Lower panel: same quantity for the 



model Eq. ( 24 \ without accounting for the hopping modula- 
tion, Ji^i+i — J =const. 



To quantitatively verify that the hopping modulation 



does not significantly affect the physics of the system, we 
have studied the first- and second-order coherence of the 



system with the full ab-initio Hamiltonian Eq. (24), and 



compared the results to the case in which we neglect the 
hopping modulation, namely Eq. (24 1 with Ji,i+i = J 



— const. The results are shown in Fig. 27 for the mo- 



mentum distribution and in Fig. 28 for the noise correla- 
tions. We observe that taking into account the hopping 
modulation does not alter the most significant qualitative 
features of these observables. As for the momentum dis- 
tribution, we observe that the hopping modulation leads 
to a stronger suppression of coherence for strong QP po- 
tentials; in the case of noise correlations, we see that the 
hopping amplitude leads to a reduction of the signal at 
the incommensurate momentum q — kmc, and even more 
significantly of the one at g = 2fcinc, with respect to the 
central peak at g = 0. For the first observation, a possi- 
ble explanation is that the hopping modulation strongly 
suppresses hopping in the minima of the potential, and 
hence the residual coherence of the particles trapped in 
those minima. As for the second observation, a reduction 
of the signal in the noise correlations is also compatible 
with a reduction of quantum fiuctuations inside the po- 
tential wells, as observed in Section |VI| when taking the 
classical limit J ^ (see Fig. 23 1. 
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VIII. SUMMARY AND CONCLUSIONS 

In summary, we have presented an extensive numerical 
study of the physics of interacting lattice bosons in a 
quasi-periodic (QP) potential, in direct connection with 
recent experiments on cold bosons in incommensurate 
optical superlattices'^. 

We have investigated the exactly solvable limit of hard- 
core repulsive bosons, where the system can be mapped 
onto n on-in teracting fermions, so that the well-known 
physicsPEslof one-dimensional particles in an incommen- 
surate superlattice directly apphes to the many-body sys- 
tem. There the emergence of gaps in the single-particle 
spectrum translates into the appearence of a gapped in- 
commensurate band insulator (IBI) phase in the many- 
body phase diagram, together with the gapless Bose glass 
(BG) phase. This picture survives also in the more real- 
istic case of softcore interactions, for which the intricate 
phase diagram is mapped out in the case of a strong QP 
potential (equalling in strength the inter-particle repul- 
sion). In particular, making use of a generalized atomic- 
limit approach, the appearence of gaps in the many-body 
spectrum is quantitatively ascribed to quantum finite- 
size effects, due to the finiteness of the potential wells in 
the QP potential. 

We have furthermore discussed the behavior of the sys- 
tem in a QP potential plus a parabolic trap, as required 
for a meaningful comparison to the experiments. The 
concept of a local phase, which is at the heart of the 
local-density approximation (LDA), applies successfully 
to a system without (pseudo-) disordered potentials ex- 
hibiting homogeneous MI and SF phases, but this con- 
cept fails to apply to the case of a BG, whose fundamental 
signature, a finite local compressibility, comes from local- 
ized gapless excitations which most likely do not appear 
in a single realization of the trapped system. A finite lo- 
cal compressibility can only be attained on average over 
the various local potential realizations (namely over dif- 
ferent spatial phase shifts) . Nonetheless the LDA picture 
for a local BG phase might be too restrictive, namely non- 
trivial low-energy particle-hole excitations might appear 
in the trap which are not present in the bulk system, 
due to accidental quasi-degeneracies of particle and hole 
states between different regions of the trap. 

The above results show that the behavior of the Bose- 
Hubbard model in a QP potential is more complex t han 
that of more widely investigated dirty-boson systemd^, 
and that in this respect a QP potential applied to the 
Bose-Hubbard model cannot be ascribed to the family 
of short-range correlated disorder potentials, given that 
correlations in such a potential are actually long-ranged. 
Both at the single-particle level and at the many-body 
level a QP potential leads to a behavior which is inter- 
mediate between that of a system in a random potential 



and that in a commensurate periodic potential: indeed 
the potential's spatial features albeit aperiodic, have a 
definite upper bound in length scale, which introduces in 
turn a characteristic energy scale in the system. Such 
energy scale shows up in the opening of gaps at incom- 
pressible regions of the phase diagram. 

If a second QP potential at a different frequency is su- 
perimposed to the first one, preliminary results show that 
the picture of a short-range correlated random potential 
is recovered, as the resulting physics of the Bose-Hubbard 
model in such a potential appears to mimic that of the 
same model in a truly random potential. In presence 
of parabolic confinement, we have moreover shown that 
a 2-color QP potential can realize a BG phase all over 
the trap on average over spatial phase fluctuations of the 
two potential components. Hence the gap between the 
single-color QP potential and an ideal (white-spectrum) 
random potential can be possibly bridged by enriching 
the QP potential with uniquely one more color. This 
conclusion has immediate practical consequences for the 
ongoing effort of realizing analog quantum simulators of 
correlated systems in a random environment within an 
optical-lattice setup, because the experimental require- 
ment to bridge the gap between pseudo-disorder and true 
disorder could be in practice the addition of one or a few 
extra laser standing waves. 

A flnal, delicate aspect is the identification of an exper- 
imental signature for the occurrence of a BG insulating 
phase. This issue appears unresolved by the currently 
accessible observables, e.g. through time-of-fiight mea- 
surements. This is intrinsically due to the fact that a 
BG is not characterized by correlations, which are gen- 
erally akin to those of more conventional insulators, but 
by the spatial and energetic structure of excitations. A 
method to resolve the nature of these excitations, based 
on measuring the system's response upon variations of 
the trapping frequency, will be presented in a forthcom- 
ing publication^. 
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